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Abstract. Thermal conductivity Kxx{T) under a field is investigated in d^2_y2-wa.ve superconductors and 
isotropic s-wave superconductors by the linear response theory, using a microscopic wave function of the 
vortex lattice states. To study the origin of the different field dependence of K^^iT) between higher and 
lower temperature regions, we analyze the spatially-resolved thermal conductivity around a vortex at each 
temperature, which is related to the spectrum of the local density of states. We also discuss the electric 
conductivity in the same formulation for a comparison. 

PACS. 74.60Ec Mixed state - 74.25Fy Transport properties - 74.25Jb Electronic structure 



, ^, 1 Introduction 



Recent advance to synthesize new exotic superconducting 
materials further requires the experimental probes to pre- 
, cisely identify their pairing functions consisting of the or- 
■ bital and spin components. The orbital part of the pairing 
' function determines the nodal structure of the energy gap 
on the Fermi surface. Thermal conductivity is one of the 
standard techniques to probe the node of the gap structure 
One can basically extract the gap topology 
such as line or point nodes by analyzing its dependence 
\ on the temperature T. 

In the vortex state under a magnetic field, the thermal 
conductivity is affected by the low energy quasiparticlc 
state around the vortex 1^,0,1). So far, the thermal con- 
ductivity in the vortex state has been investigated by the 
\ theory of gapless superconductors at high field P, pT[ | , o r 
the theory considering vortices as scattering centers pT| . 
Recently, the thermal conductivity in the vortex state of 
high-Tc superconductors attracts much attention, because 
' the quasiparticlc states are qualitatively different in the 
vortex state between the d-wave pairing case and the con- 
ventional s-wave pairing cases. In the s-wave pairing, low 
energy quasiparticle states are bounded around the vor- 
tex core p^ , |l3| . In the d-wave pairing, low energy quasi- 
particle states around the vortex extend outside of the 
vortex core due to the line node of the superconducting 
gap @|l|,|l|,|l|,[8[|l|l. The contribution to the thermal 
conductivity from the quasiparticles outside of the core 
is investigated by the theory of the Doppler shift (or the 
Dirac fermion) in the d-wave pairing case |2^,^,^. Since 
these theory neglect the contribution from the quasipar- 



ticle within the vortex core they are not applied to 
the s-wave superconductors. 

Here, we calculate the T-dependence of the thermal 
conductivity under a field. Our calculation is based on a 
microscopic theory of Bogoliubov-de Gennes (BdG) equa- 
tion for describing the vortex state |^,^|2^ and stan- 
dard linear response theory for the thermal conductiv- 
ity [ p6| , |2^ , assuming a clean-limit type II superconduc- 
tors. CJur calculation includes all contributions from the 
inside and the outside region of the vortex core. The spa- 
tial distribution of the thermal conductivity is calculated 
from the wave function of the vortex lattice state, and an- 
alyzed by compared with the spatial distribution of the 
local density of states (LDOS). 

The rest of this paper is organized as follows. In Sec. 
H, we describe our formulation based on the BdG equation 
and the linear response theory. In Sec. ||, we study the de- 
pendence of the thermal conductivity on the temperature 
in the vortex state. We also show the position-resolved 
thermal conductivity and its energy decomposition, and 
discuss the relation with the LDOS. In Sec. ^, we investi- 
gate the electric conductivity in the same formulation, and 
discuss the difference of the quasiparticle contribution be- 
tween the thermal conductivity and electric conductivity. 
The summary and discussions are given in Sec. I 



2 Formulation 

2.1 Bogoliubov-de Gennes equation 

We obtain the wave function in the vortex lattice state 
by solving the BdG equation for the extended Hubbard 
model [E3,24,23] defined on a two dimensional square lat- 
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quasi particle structure as previous theoretical studies both 
for the s-wave [^,13 case and for the d-wave ||l^,|l5|,p^ 
|l^,|l8|,|l9) case. Here we briefly discuss the BdG equation 
for the s-wave pairing and the d-wave pairing cases, which 
is written as 



I D^^ -K* 



— Ea 



where 



Ilj = Uj expp- 



A(r).dr], 



(1) 

(2) 
(3) 

(4) 



with the on-site interaction U, the flux quantum (f>o and 
the chemical potential fi. The nearest neighbor (NN) trans- 



fer integral t 



t and the NN interaction Vi 



'1,3 



V 

for the NN site pair and r,i±e- The vector potential 
A(r) = iH X r in the symmetric gauge. Since we assume 
an extreme type-II superconductor, the internal field term 
of A(r) is neglected. The self-consistent condition for the 
pair potential is 



A. 



2 ^ 



u„(r,)<(rj)tanh(£;„/2T). 



The periodic boundary condition is given by the sym- 
metry for the translation R = ^^R" + ^j,R.y, where R° = 
(aAf,.,0) and R° = (0,aiVr) are unit vectors of the unit 
cell for our calculation. Then, the translational relation 
is given by ^^(r -I- R) = UQ(r)e*^('''^)/^, ■DQ(r -I- R) = 
w„(r)e-^x(""^^)/^ Here, 



O-jr Ojr 

x(r,R) = -— A(R)-r-27rZ,(/,-?j,) + — (Hxro)-R (9) 

in the symmetric gauge. The vortex center is located at 
ro -I- |(3R° -t- R°). The phase factor in Eq. (|) is nec- 
essary to satisfy the translational relation Ad{v + R) = 
Z\<j(r)e'x('-,R). 

The following parameter values are chosen in the cal- 
culation. The average electron density per site ~ 0.9 by 
appropriately adjusting the chemical potential ^. We nor- 
malize all the energy scales by the transfer integral t. For 
the s-wave case, we set U — — 2.32t and V = 0. The re- 
sulting order parameter Aq — 0.5t at T = = 0, and 
the superconducting transition temperature Tc ^ 0.27t. 
For the d-wave case, we set J7 = and V — — 4.2i. Then, 
Ao = l.Ot, and Tc ^ QA2t. Our results do not qualitatively 
depend on the choice of these parameters. 



(5) 2.2 Local density of states 



The s-wave pair potential is given by 
A,{vi) = UA,^,. 
The (i2.2_j,2-wave pair potential is 
V 



A. 



y,i- 



(6) 



^-S,.) (7) 



with 



A±e^i = Z\,,,;±e exp[i — 
00 



(ri+ri±,)/2 



A(r)-dr]. (8) 



In order to calculate physical quantities, we must con- 
struct the Green's functions from Ea,Ua{v), Vai'r) in the 
formulation of imaginary time r and Fermionic Matsub- 
ara frequency a;„ = 2nT{n + i). The Matsubara Green's 
functions is given by 



g{r, r',iw„) 



_ f gi2(r,r',i 



(10) 



g2iir,r', iLOn) .g22(r, r', i^n) 
where matrix components are Fourier transformation of 



We study the case of the square vortex lattice where 
the NN vortex is located in the direction of 45° from 
the a-axis. This vortex lattice configuration is suggested 
for d-wave superconductors and s-wave superconductors 
with fourfold symmetric Fermi surface |l|,||,|9[|^,|l|. 
The unit cell in our calculation is the square area of iV^^ 
sites where two vortices are accommodated. Then, H = 
2(j)o / {aNr)'^ with the lattice constant a. Thus, we denote 
the field strength by as H^^ . Since H should be com- 
mensurate with the atomic lattice, our formulation does 
not treat the field dependence continuously. We consider 
the area of iV^^ unit cells. By introducing the quasi-momentum 
of the magnetic Bloch state, k {2Tr /aNrNk){lx,ly) ■ 
{lx,ly = 1,2,- • •, A^fc), we set Ua{r) = Ua{r)&^'' ,Va{r) = 
Va{j^)e^^''^ ■ Then, the eigen-state of a is labeled by k and 
the eigenvalues obtained by this calculation within a unit 
cell. 



gn(r,r,r',r') = -{TA^j^ir,T)4>\ir' ,t')]), 
gi2(r,r,r',r') = -(TAi^^{v,T),p^{r' ,t')]) , 
g2i(r,r,r',r') = -(r^iV-jCr, r)V| (r', r')]), 
522(r, T, r', r') - -{TAi^lir, r)V'i(r', r')]). 

The Green's functions in Eq. ( p^ ) as follows, 

"a(rX(r') 



gii(r,r',itj„) 



5i2(r,r',iu;„) = 



E 

a 

E 



ILOn — Ea 

Ma(rX(rO 
iujn — E„ 



521 r,r ,icj„ = > - ^ ' 



522(r,r',ia;„) = ^ 



lUJn — Ea 
Va{Y)v*a{Y') 
ILUn — Ea 



(11) 

(12) 
(13) 
(14) 



(15) 
(16) 
(17) 
(18) 
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From the LDOS is given by the thermal Green's func- 
tions as 



[E, r) = - -Imgn (r, r, icj„ ^ S + iO^ 

TT 

= Y.\u^{v)\H{E - E^) 



for the up-spin electron contributions, and 

Ni {E, r) = -lmg22 (r , r, -icj„ E + iO^ 



(19) 



with a — (a,0). Four electron field operators in Eq. (|2£ 
are decomposed by using the Matsubara Green's functions 
in Eqs. (Ill)-® . 

For the study of the thermal conductivity, we have to 
introduce the dissipation term 77 in the Green's function. 
Then, the retarded and advanced Green's functions are, 
respectively, given by G^(r, r' ,lu) = g{r, r', icj„ oj + irj) 
and G-^{r,r',uj) = ^(r, r',iu;„ — > uj — ii]). Therefore, in the 
spectral representation, we obtain 



J2Kir)\'SiE + E^) 



(20) 



5(r,r',iw„) = 



duj A{r, r', uj) 
27r iuJn — w ' 



(26) 



where 



for the down-spin electron contributions. Therefor, the 
LDOS is given by 

NiE,r) ^N^{E,r)+N^{E,r) 

= ^{\u^ir)\HiE - E^) + \v^{v)\H{E + E^)}. 



i(r, r', c.) = (Gi^(r, r', c.) - ^^(r, r', 

Ma(rX(r') Wa(rX(r') 
(r)<(r') 



(27) 



(21) 



with 5r,{oj) = 77[7r(w2 r;^)]-!. 

As a result, the heat-heat correlation function is re- 
duced to 



2.3 Linear response theory 



We calculate the thermal conductivity following the method 
of Refs. |2^ and |^ . According to the linear response the- 
ory, thermal current hxiji) flowing to the a;-direction at 
ri-site is given by 

^ ^Re{ j^(5a;a:(ri,r2,iJ7„ ~> Q + iO+)}r2-.o 



d!2 



--^'^F{Ea,Ep) 



Am? 



x(-V,T(r2)), 



(22) 



X [PiUQ(ri)M^(ri') + Uairi)P\,ul{rv) 
-HPlw„(ri)w^(riO + i;„(ri)Pr'y^(riO] 
x[P2UQ(r2)M^(r2') + MQ(r2)P^,u^(r2') 
+Pjz;„(r2)z;;(r20 + t^c(r2)P2't';(r2')]|i=i',2=2', (28) 



when the small temperature gradient — Va;r(r2) along the 
a;-direction is applied at r2-site. The heat-heat correlation 
function is deflned by 



where 



Qxx{riTi,r2T2) = {Tr[hx{riTi),hx{r2T2)]) 



-if2„(ri-T2; 



Qxxiri,r2,if2n) 



(23) 



F{Ea,Ep) 

dw r dw' 
2t: J ^ 

,a;V(c^) - ^'V(^') 



X P- 



S,^{UJ - Ea)S,j{Lj' - Efj) 

iTruj^f'{Lu)S{oj-Lu')]. (29) 



in the formulation of imaginary time r and Matsubara fre- 
quency fin- The heat current operator h(r^-,r) is written 
as 



When the temperature gradient ViT(r2) is uniform, 
the position-resolved thermal conductivity is written as 



Kxx (1*1 ) 



hx{ri) 



V4 (r,' , r' ) V. (rj , r) (24) 



in terms of the electron field operators '0cr(rj,T). The 
a;-component of the momentum operator Pj in the dis- 
cretized square lattice is defined as 



[P,Mr„r)]x = a^e'(-/^'')-^C-^+-/2)V'.(r, + a, r) 

(25) 



- ^M^^ ^ Q..(ri, r2, inn^f2 + i0+)}.(30) 
The spatially averaged thermal conductivity 

Kxx = J^'^Kxxiri) (31) 
ri 

is observed in the experiment. N is the total number of 
sites. At zero field, our formulation for Kxx is reduced to 



4 



Please give a shorter version with: \authorrunning and \titlerunning prior to \maketitle 




the well-known formula for the uniform superconductors 
in the presence of impurity scattering p6| . 

Using the wave function Ua{r), Wa(r) and the eigen- 
energy Ea in Eq. we calculate the dependence of Kxx 
on the temperature and the field. And to analyze this 
behavior, we also study the spatial structure of k,j;x{y), 
i.e., the local contribution to Kxx, Here, we neglect the 
principal value integral term KeF^Ea, Ep) in Eq. (p9|), 
because the contribution from this term vanishes in the 
spatial average of Kxx{j^)- We typically choose rj ~ O.Olt. 



3 Thermal Conductivity 

3.1 Temperature dependence of thermal conductivity 

The T-dependence of Kxx{T) is shown in Figs. |l](a) and (b) 
for the s-wave and the d-wave pairing cases, respectively. 
For the s-wave pairing, it shows exponential T-dependence 
due to the full gap of the s-wave superconductivity at H — 
0. It changes into a T-linear behavior at low T region in the 
vortex state at H ^ 0, reflecting low energy quasiparticle 
state around the vortex. We see also that the deviation 
from the expected T dependence occurs by the impurity 
effect at very low temperature around T ^ 2r] = 0.02t. 
As for the d-wave pairing in Fig. |l|(b). The zero field case 
shows the T^-behavior , which ischaracteristic of the 
line node of the dx^-y^-'wave superconductivity. For H ^ 
0, it is modified to the T-linear dependence at low T 
region. 

It is seen in Fig. |l| that there exists a crossover temper- 
ature T* both in the s-wave and the d-wave pairing cases, 
when we consider the dependence of k,xx{T) on the mag- 
netic field. At lower temperature T < T* , Kxx{T) increases 
with raising magnetic field. However, at higher tempera- 
ture T > T*, Kxx{T) decreases as a function of H. It is 
also noteworthy that Kxx{T) shows the T-linear behavior 
at T < T*, while it deviates from T-linear at T > T*. 
In our parameter, T* ^ O.lOt in the s-wave case, and 
T* ^ 0.06i in the d-wave case. 



3.2 Position-resolved thermal conductivity 

To understand the abovementioned difference between T < 
T* and T > T*, we analyze the position-resolved Kxx{j^) 
of Eq.(pO|), and investigate how the local thermal fiow 
contributes to the total thermal conductivity. The spa- 
tial structures of K.xx{r) are shown in Fig. ^ for the s-wave 
pairing and in Fig. ^ for the d-wave pairing from low tem- 
perature to high temperature. As for the vortex core size, 
the order parameter recovers at two or three lattice sites 
from the vortex center at low temperatures. It is seen for 
both pairing cases that the heat flows exclusively at the 
core region at low T. It comes from the available low en- 
ergy excitations around the vortex core. Around the core 
region the heat flow extends to the NN vortex direction. 
It is due to the inter-vortex quasi-particle transfer effect. 
With increasing T, the the contributing region of Ka;x(r) 
becomes wider around the vortex core and the lines be- 
tween NN vortices, as seen in Figs, ^(b) and ||(b). On the 
other hand, at higher temperature T > T*, the dominant 
contribution comes from the outside region of the vortex 
core. Since Kxx(j) is suppressed at the vortex core, the 
vortex core behaves as if it is a scattering center for the 
thermal flow. Also along the line connecting NN vortices, 
i-^xxi^) is slightly suppressed. At further high T, the heat 
flows rather uniformly at all region. These higher temper- 
ature structure of Kxxi'<^) comes from the contribution of 
the scattering state at £^ > |Z\|, as discussed later. 

Let us discuss important differences between the s- 
wave and the d-wave pairing cases. At lower temperature 
T — O.Oli, Kxx{j^) is well localized at the core region in 
the s-wave case (Fig. ||(a)) compared with the d-wave case 
(Fig. ^(a)). This difference comes from the line node con- 
tribution of the d-wave superconductivity. Due to the line 
node, the low energy quasiparticle states extend outside of 
the vortex core, especially to the 45° direction from a and 
b axis of the crystal lattice At higher temperature 

T = 0.09t, Kxx{j^) is well suppressed at the core region in 
d-wave case (Fig. ||(c)), and it is not too suppressed at the 
core region in s-wave case (Fig. |3|(c)). To understand these 
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(a) 



(b) 



(d) 



Fig. 2. The vector plots of the nxxi^) in the s-wave pairing case at T/t=0.02 (a), 0.05 (b), 0.09 (c) and 0.13 (d). H 
The vortices are located at the center and the four corners in the figure. Arrow size is proportional to flow strength. 
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(b) 




(d) 



Fig. 3. The vector plots of the ^^^^^(r) in the d-wave pairing case at r/t=0.02 (a), 0.05 (b), 0.09 (c) and 0.13 (d). H 
The vortices are located at the center and the four corners in the figure. Arrow size is proportional to flow strength. 




Fig. 4. Ea- and i5^-dependence of the function \\aF(Ea, Ep) 
of Eq. (H). In this figure, we set r;/t=0.01 and T/f=0.10. 

characteristic behavior of n^x (r) , we consider the relation 
between the thermal conductivity and the LDOS. 

3.3 Relation with the local density of states 

In Eq. ( psi) , temperature dependence comes from the func- 
tion F{Ea.,Ep) in Eq. (|29|). It determines which energy 
level dominantly contributes to the thermal conductiv- 
ity. To see it, we show the Ea and Efs dependence of 
lviiF{Ea,Ep) in Fig. ^. From the figure, we see that ImF(£', 
becomes large on the line Ea — E/j. Along the line Ea — 
Efs = E, we obtain 

lmF{E,E) ^ E^f{E). (32) 



Then, ImF(£', E) vanishes at i? = 0, and has two peaks at 
finite energy, which are symmetric with respect to = 0. 
The distance of these peaks is about bV^T. It means that 
the quasiparticle states at these peak energy E ^ ±2.5T 
dominantly contributes to Kxx(j)- At low temperature, the 
dominant contribution comes from the low energy quasi- 
particle state around the vortex core. However, at higher 
temperature, the scattering states at i? > A{T) domi- 
nantly contribute to the thermal conductivity. In this re- 
spect, thermal conductivity is qualitatively different from 
other physical quantities such as electric conductivity, spe- 
cific heat [|l|,Q|2|], nuclear magnetic relaxation time [E4l 
psf . In these quantities, the quasiparticles at i? ~ gives 
largest contribution in all temperature regions. We discuss 
the electric conductivity at Sec. ^. 

At each energy level Ea , the contribution to the spatial 
structure of Kxx{^) is determined from the spatial distri- 
bution of the wave functions Uai^) and Wa(r). Then, we 
show the spatial structure of the LDOS at several energies 
for the s-wave pairing in Fig. ^ and for the d-wave pairing 
m Fig. I In the s-wave pairing, the low energy quasiparti- 
cle states are bounded within the vortex core region, and 
the quantized energylevels appear at half integer energy 
En = {n + ^)E/<i Here, Ea is the level spacing 

of the order Aq/Ep. Aq is the superconducting gap at 
zero field and Ep is Fermi energy. With the BdG theory 
^ ^lean limit, since the vortex core radius shrinks to the 
atomic scale with lowering temperature by the Kramer- 
Pesch effect , the quantization effect eminently ap- 

pears at T ~ 0. Then, there are no states just at i? = 
because of the small gap by the quantization in the s-wave 
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pairing. At the lowest energy level - 0.17t in Fig. |(a), 
the LDOS N{E, r) has sharp peak at the vortex center. It 
is a bound state in the vortex core. At higher energy, the 
LDOS has peak along a circle around each vortex (Fig. 

(b)). The radius of the circle increases with raising en- 
ergy. We also see the small ridge between NN vortices. It 
is due to the inter- vortex transfer of the low energy bound 
states. At E ^ A, the circle of the LDOS peak around the 
vortex overlaps each other (Fig. ||(c)). At higher energy 
than A, the LDOS are reduced to the uniform structure, 
though the LDOS is slightly suppressed at the vortex core 
region (Fig. 1(d)). These energy dependence is consistent 
with the results of the quasiclassical calculation p^,p5[. 

For the rf-wave pairing in Fig. ^ since the low energy 
states extends outside of the vortex core due to the node 
of the dj,2_j,2-wave superconducting gap, energy levels be- 
comes continuous [|8,|l^,^,|^|2^. In Fig. Q(a), there is 
a peak of the LDOS at the vortex core aX E = 0, which 
corresponds to the zero-bias peak in th spectrum at the 
vortex center. With increasing E, the peak of the LDOS 
is shifted to the outside of the vortex, and it is slightly 
suppressed at the vortex center, as shown in Fig. ^(b). In 
the d-wave case, the large LDOS region shows the four- 
fold symmetric structure instead of the circle structure of 
the s-wave case. At E A, the LDOS peak around the 
vortex is shifted to the boundary region between vortices 
(Fig. ^(c)). At higher energy than A the LDOS are re- 
duced to the uniform structure (Fig. ||(d)). 



3.4 Energy decomposition of thermal conductivity 

To discuss the contribution of the LDOS structure to the 
spatial structure of Ka;a;(r), We decompose Kxx{^) of Eqs. 



and (|30| ) into the low energy contribution from \ Ea\, \ Ef}\ 
A{T) /2 and the high energy contribution from \Ea \ jjEp \ > 
A{T)/2. The s-wave pairing case is shown in Fig. The 
upper panels (a) and (b) present the low energy contri- 
bution of Kxxi''^), which are localized around the vortex 
core and along the line connecting NN vortices. It is be- 
cause the low energy states for \E\ < A{T) are localized 
around the vortex core and there are some inter-vortex 
quasiparticle transfer along the NN vortices direction, as 
shown in Fig. ||(a) and (b). The lower panels (c) and 
(d) show the higher energy contribution of Kxxi''^) from 
I^^qI, \Ep\ > A{T)/2. In this energy range, the wave func- 
tions are dominantly located outside of the vortex core, as 
shown in Fig. |^(c) and (d). Then, the higher energy con- 
tribution of Ka;a;(r) is suppressed around the vortex core. 
The suppression along the NN vortices directions shown 
in Fig. Wd) also reflects the spatial structure of the LDOS 
in Fig. g(d). At low temperature in Fig. 0(c), spatial struc- 
ture is determined by the wave function at E' ~ A{T)/2 
in Fig. ^(b). However, at high temperature case in Fig. 
0(d), the contribution of the wave function at E > A{T) 
in Fig. ^(d) is dominant. 

The d-wave pairing case is shown in Fig. |^. In the up- 
per panels (a) and (b) for the lower energy contribution, 
'^2;x(r) is larger around the core and the lines connecting 
NN vortices. It is broadly extending around the vortex, 
compared with the s-wave case, because the wave func- 
tions are also broadly extending around the vortex core 
as shown in Fig. ^(a). At higher temperature, Kxx{'>^) is 
large outside of the core in Fig. ||(b), though LDOS at 
E — A{T)/2 are little localized around the core. In the 
lower panels (c) and (d) for the higher energy contribution, 
'ta:£c(r) is large outside of the core at every temperature. 



< 
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Next, we investigate the weight of the energy-decomposed 
contribution for the spatiaUy averaged k^^. The temper- 
ature dependence is presented in Fig. ^. For both pair- 
ing cases, we can see that the low (high) energy contri- 
bution is dominant at low (high) temperature. The low 
(high) energy contribution of the d-wave pairing case is 
larger (smaller) compared with the s-wave pairing case. 
It is because the d-wave pairing case has larger DOS at 
\E\ < A{T) because the low energy excitation widely ex- 
tends outside of the vortex core region due to the line 
node of the superconducting gap. These DOS difference 
between the d-wave and the s-wave pairing cases is also 
shown by the quasiclassical calculation (Fig. 18 of Ref. 
IIitI ) . The very low temperature behavior in Fig. ^ (a) in 
the s-wave pairing at T < 10~^ reflects the small gap of 
the quantized energy level in the s-wave pairing. 



4 Electric conductivity 

We can also calculate the electric conductivity, if we con- 
sider the electric current operator instead of the thermal 
current operator. Following the same procedure in Sec. 
2.3, the position-resolved electric conductivity is given by 




Fig. 7. Energy decomposed Kxx{r) in the s-wave case at T/t — 
0.02 [(a),(c)] and 0.09 [(b), (d)]. (a) and (b) are for the low 
energy contribution from \Ea\, |-E/3| < A{T)/2. (c) and (d) are 
for the high energy contribution from \Ea\, j-B/al > A(T)/2. 



d 1 



a,,(ri) =Im{ — -^QSUri,r2,il?„ ^ 12 + iO+)}|,,^o 



dQN 

with the correlation function 

Qtir,T,,V2T2) = (r,[j,(riTi),J,(r2T2)]) 

= r^e-'^^"(-^--^)QSUri,r2,iX2„) 

n 

of the electric current operator 



(33) 



(34) 



(35) 



Using the wave functions of the BdG equation, we obtain 




af3 



x[PiUairi)ul{ri,) + Uairi)P\,Up{rr) 
-Pl^;„(ri)i;;(riO - t;„(ri)Pi,«;(riO] 

X [P2Ma (r2 )W^ (r2' ) + M„ (1-2 )P2' Up (r2' ) 

-P^iv(r2)t;;(r2') - «a(r2)P2't;;(r2')]|i=i',2=2', (36) 



where 



du) f doj' 



2-n J 2tt 



Sniuj - Ea)Sr,{uj' - Efi) 



X P^^^Mfe^ + i^/'(^)'5(^ - u;')]. (37) 



Fig. 8. Energy decomposed n^xijc) in the d-wave pairing case, 
at T/t = 0.02 [(a),(c)] and 0.09 [(b), (d)]. (a) and (b) are for the 
low energy contribution from |-E,3| < A(T)/2. (c) and (d) 
are for the high energy contribution from |-Ea|, \Ei) \ > A{T)/2. 



We show the Ea and Ef3 dependence of Inii^*^' (£'„ , Ejs) in 
Fig. |o[ The principal value part of Eq. (^^ vanishes by 
the spatial average of axx(ji). Then, we consider only the 
contribution ofli[nF'^^{Ea,Ep), neglecting ReF'^^{Ea, E/3). 
The function lmF°^{Ea, E/^) becomes large for Ea ^ Ep. 
Along the line E^ ^ Ep ^ E, we obtain ImF''\E,E) = 
f'{E). It has maximum at i? = in all temperature re- 
gion. Then, the low energy quasiparticle states around the 
vortex core dominantly contribute to a^x (r) even at higher 
T. We present (Ta;a;(r)_in Fig. 0(a),(b) for the s-wave pair- 
ing case and in Fig. ll[c),(d)for the d-wave pairing case. 
Even at higher temperature [(b) and (d)], the dominant 
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H=H2o 

\E^p\<A(jy2 

\E„,p\>A(J)/2^ 



T/t 



T/t 



Fig. 9. Temperature dependence of the energy-decomposed Hxx- The low energy contributions from \Ea\, < A{T)/2 and 
the high energy contributions from \Ea\, > A{T)/2 are shown with the total k^i! (Solid line). H = H2o- (a) s-wave case, 
(b) d-wave case. The arrows Ta-Td show T/t— 0.02, 0.05, 0.09 and 0.13, respectively. We show Kxx{r) in Figures |2|, ^, |^ and ^ at 
these temperatures. 



Im F^'(E„,E^) 




Fig. 10. Ea- and i5/3-dependence of the function 
lmF'''{Ea, Ep) of Eq. @. In this figure, we set r)/t=0.01 and 
T/t=0.10. 



contribution to a^x conies from the vortex core region as 
in the low temperature case [(a) and (c)]. This is con- 
trasted with the thermal conductivity case, whose dom- 
inant contribution comes from the outside of the vortex 
core at high temperatuer, as discussed in Sec. |[ Compared 
to the s-wave pairing case[(a) and (b)], Uxx{t^) widely ex- 
tends toward the outside of the vortex core in the d-wave 
pairing case [(c) and (d)]. It reflects the extending low 
energy quasiparticles due to the line node of the d-wave 
superconducting gap. 



5 Summary and Discussions 

We have formulated thermal conductivity in mixed state 
based on a microscopic theory of BdG equation and linear 
response theory. The T-dependence of thermal conductiv- 
ity K,xx for the s-wave and the d-wave pairing is calculated. 
Their behaviors are analyzed in terms of the position- 
resolved thermal conductivity Hixx{^)- And we discuss the 
relation between Kxx{Ti^) and the LDOS of the quasiparti- 
cles around the vortex. 
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(c) (d) 

Fig. 11. The vector plots of the position-resolved electric 
conductivity axx{r) in the s-wave pairing case at T/t=0.02 (a) 
and 0.09 (b), and in the d-wave pairing case at T/t— 0.02 (c) 
and 0.09 (d). H = H^o- The vortices are located at the center 
and the four corners in the figure. Arrow size is proportional 
to flow strength. 



There is a crossover temperature T*. At lower tem- 
perature T < T* , Kxx is increased with raising magnetic 
field. In these temperature region, thermal flow is domi- 
nantly carried by the low energy quasiparticles around the 
vortex core and their inter- vortex transfer. Then, Ka;a;(r) 
is large around the vortex core and the lines connecting 
NN vortices. On the other hand, at higher temperature 
T > T*, Kxx is decreased at higher field. In these temper- 
ature, the contribution to the thermal conductivity comes 
from higher energy quasiparticles including the scattering 
state at i? > A{T). Therefore, Kxx{ii^) is suppressed at 
the vortex core region. Then, vortex core works as if the 
scattering center for the thermal flow. These contributions 
from the higher energy states at higher temperature is a 
characteristic of thermal conductivity. For other quantities 
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such as electric conductivity, specific heat, nuclear mag- 
netic relaxation time, the low energy state gives largest 
contribution at all temperature regions. 

The difference between the s-wave pairing and the d- 
wave pairing comes from the node structure of the super- 
conducting gap. The d-wave pairing case has larger contri- 
bution from the low energy quasiparticle states. The low 
temperature distribution of K^xi''^) is broadly extending 
around the vortex in the d-wave case, because the wave 
functions are broadly extending around the vortex core 
due to the node structure. 

Our caluclation has reproduced experimental results 
of the T-linear behavior at low temperature in the vortex 
states {H ^ 0), and the existance of the crossover tem- 
perature T* in the field dependence. The crossover tem- 
perature T* are reported both in the conventional s-wave 
superconductor such as Nb (Ref. 0) and in the high-Tc 
superconductors such as Bi2Sr2CaCu208 (Ref. and 
in the /-wave superconductor UPts (Ref. Q). These ex- 
perimental results have been explained as follows. At low 
temperature, vortex assist the thermal flow due to the 
low energy quasiparticle state around the vortex core. At 
higher temperature, vortex behaves as the scattering cen- 
ter for the thermal flow. Our numerical results of n^xij) 
are consistent to this picture. While the field dependence 
of n,xx is important, our calculation cannot examine the 
contunuous field dependence because we consider the field 
H = 2(t)Q / {aN r)^ depending on the size of the unit cell. 
At higher temperature, there appears the effect of the T- 
depending r] due to the inelastic scattering by antiferro 
magnetic spin flactuations |^,^|3^. For the further ex- 
tention of our calculation, we will examine the effect of 
the T-dependence or the position dependence (inside or 
outside of the vortex core) of the scattering parameter rj. 
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